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We study theoretically the formation spectra of deeply bound pionic atoms expected to 

y I be observed by experiments with high energy resolution at RIBF/RIKEN, and we discuss 

in detail the possibilities to extract new information on the pion properties at finite density 

from the observed spectra, which may provide information on partial restoration of chiral 

symmetry in medium. We find that the non-yrast pionic states such as 2s are expected to be 

seen in the (d,^He) spectra, which will be helpful to reduce uncertainties of the theoretical 

calculations in the neutron wave functions in nucleus. The observation of the 2s state with 

QQ ' the ground Is state is also helpful to reduce the experimental uncertainties associated in 

the calibration of the absolute excitation energy. We find that the nuclear densities probed 

Qs , by atomic pions are quite stable and almost constant for various atomic states and various 

1^^ ■ nuclei. Effects of the pion wave function renormalization to the formation spectra are also 

• ' evaluated. 

!>■ 

^—^ ■ §1. Introduction 

Properties of hadrons at finite density and temperature are extremely interesting 
in the contemporary hadron-nuclear physics because they provide important hints 
^ , to explore the relation between the symmetry breaking pattern of QCD and the 

C^ ' observed hadron properties for the study of the QCD vacuum structure. 'I'' One of 

the good systems to observe in-medium hadron properties is the deeply bound pionic 
atomP ™ which is a 7r~ atomic bound state hardly observed by X-ray spectroscopy, 
such as the Is or 2p states in heavy nuclei. The deeply bound states have been 
experimentally produced in the ((i,^He) reactions with Pb and Sn isotope targets 
at GSI^T^ by following theoretical predictions .'^"'1^ In the latest experiment,'^ 
the energy shifts and widths of the Is states have been precisely measured in three 
Sn isotopes and isospin-density dependence of the s-wave pion-nucleus potential has 

typeset using VTpT^.ds (Ver.0.9) 
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been deduced. From these observations, reduction of the chiral order parameter {qq) 
in nucleus was concluded. A recent model independent theoretical analysis supported 
the way to extract the in-medium quark condensate from the pionic atom data 
and showed a relation connecting the in-medium quark condensate to the hadronic 
observables.'^ For further studies of in-medium pion properties, new experiments 
were proposed to make high precision spectroscopy of pionic atoms systematically 
in RIBF/RIKEN.1^'1^ Thus, the pionic atom can be one of the best systems to 
deduce the quantitative results for the meson properties and the partial restoration 
of chiral symmetry in medium around normal nuclear density at zero temperature. 
Based on these theoretical and experimental developments, we think that we 
should consider now the possible future directions of the studies of the deeply bound 
pionic atoms after 15 years from the discovery of the deeply bound pionic atoms at 
QgjpJ,[5J,[T6Ji,[T7|i \Yg discuss the following points in this article in detail, 
(i) Predictions of the pionic atom formation spectra by the ((i,^He) reactions on 

Sn and Te isotope targets proposed in Refs. |Ilj),[l5j) at RIBF/RIKEN, 
(ii) Possibility to observe pion properties and to determine the value of the chiral 

condensate at different nuclear density from p = 0.6/9o to obtain information on 

the symmetry breaking parameters beyond the linear density approximation, 
(iii) Possibility to determine the wave function renormalization factor introduced in 

Refs. fT3l) and [T8]l from pionic atom observables, 
(iv) Uncertainties included in theoretical calculations used to evaluate formation 

cross sections, 
with paying the attention to the advantages of the simultaneous observation of the 
non-yrast 2s bound state together with the deepest Is state for the same nucleus in 
the new experiments with better energy resolution.'^''^ 

In Section 2, we summarize the theoretical formalism used to connect the pion 
properties to the order parameter of the chiral symmetry and the wave function 
renormalization factor. We also mention the sensitivities of atomic pion to nuclear 
densities and the strong correlation between potential parameters. In Section 3, we 
show the numerical results and discussions for the structure and formation of deeply 
bound pionic atoms. Section 4 is devoted to the conclusion. We summarize the 
supplementary formalisms and numerical results in Appendixes. 

§2. Formalism 

We introduce briefly the formula reported in Ref. 13) in this section as the guide 
to deduce the information on chiral symmetry from the pionic atoms observables. 
We also mention the sensitivities of atomic pions to nuclear density and the strong 
correlation of the potential parameters. The theoretical formula used in this article 
to calculate the structure and the formation cross sections of the pionic atoms are 
summarized in Appendix Rl 

2.1. Chiral Dynamics of Pionic Atoms 

The spontaneously broken chiral symmetry in vacuum is expected to be restored 
partially in nuclear medium. The partial restoration of chiral symmetry takes place 
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with effective reduction of tlie chiral quark condensate {qq) in medium.'^ 'l^ Ex- 
perimental observation of the reduction of the quark condensate is not so trivial, 
since the quark condensate is not a direct observable. In Ref. [Tj) K. Suzuki et al. 
made use of in-medium extensions of the two well-known relations, Gell-Mann-Oakes- 
Renner (GOR) relatioiP^ and Tomozawa- Weinberg relation (TW),'^''^ to connect 
the quark condensate and the observables in pionic atoms, after having extracted the 
values of the parameters in the optical potential from the observed energy shifts and 
widths of the deeply bound atomic states. The 61 parameter of the optical potential 
in Eq. ()A-5|1 corresponds to the effective scattering length between pion and nucleus, 
which can be expressed by an in-medium extended TW relation: 



ry = -4ne,h = —,, (2-1) 



where /* is an in-medium pion decay constant. Assuming that the pion mass does 
not change in nucleus, one could have an in-medium GOR relation 

miff = -2m,{qq)p, (2-2) 

where ruq is the isospin-averaged quark mass mg = {niu + mrf)/2 and {qq)n is the 
in-medium quark condensate. Using these two relations, K. Suzuki et aly^ could 
get a connection between the experimental observation and the in-medium quark 
condensate through bi and /*. 

The experimental observations of the deeply bound pionic states and the at- 
tempt to connect the observables to {qq)p based on the simple extension of the 
in-vacuum relations have stimulated theoretical works to give stronger foundations 
of the analysis. It was shown based on chiral perturbation theory^-^ and on a cor- 
relation function analysis™ that the in-medium TW relation Eq. (|2-ip can be valid 
in the linear approximation of the isovector density, which leads to 

Lfree / £t 

1 - / -^^ 1 (2.3) 



^1 V/vr 

with the in-vacuum isovector vrA'^ scattering length b^^, the in-vacuum pion decay 
constant /^r and the time component of the in-medium pion decay constant /* . Fur- 
ther, it was found in a model-independent argument based on the operator relation 
that there is a sum rule for the in-medium quark condensate.'ISJ' This sum rule can 
be simplified at low density limit and gives a new scaling relation 

Zf^^ ( ¥] > (2-4) 



im) " \f: 

where Z* is the wave function renormalization for the in-medium pion (see also 
Ref. [M|) ). The density dependence of the wave function renormalization can be 
estimated at low density limit through the vrA^ scattering amplitudeI21' : 

Z: = l-I3p, (2-5) 
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with /3 = 2.17 ±0.04 fm^. Combining Eqs. ()2-3p and (|2-4p . one obtains a connection 
between the in-medium quark condensate and the experimental observables as 

The hi parameter in Eq. ()2-6p is obtained by the pionic atom data, and the rest in the 
right hand side, h^l'^'^ and Z*, are evaluated by the vrA^ scattering. For completeness, 
it is very good if one can determine also the wave function renormalization from 
the deeply bound pionic atoms instead of using the in-vacuum vrA^ scattering. We 
discuss the possibility to determine the pion wave function renormalization factor 
from the observation of the pionic atom states in Section 3.3. We also mention here 
that the estimation of the higher order effects of the density from experimental data 
is also interesting and important to explore the behavior of {qq)p beyond the linear 
density approximation.'^T^ The discussion of the possibility to deduce the higher 
order effects from the observation will be given in Section 13.21 



2.2. Nuclear densities probed by atomic pion and Seki-Masutani correlations 

As we will see in detail in Section \'6.2\ the pion in atomic states observed by the 
X-ray spectroscopy are known to be only sensitive to narrow range of nuclear den- 
sity, which is almost independent of the nuclides and the atomic states.l^ Hence, 
the structure of the pionic atoms is essentially determined by the optical potential 
strength at the effective nuclear density pe probed by atomic pion. Namely, the 
series of optical potentials which have the same potential strength at p = pe equiva- 
lently provide the almost same structure of the pionic atoms. As a consequence, we 
have a certain relation between potential parameters of the optical potentials which 
reproduce the atomic data well. We consider the s-wave part of the optical potential 
Vs as an example which includes the important piece 6i to deduce the {qq)p value 
and plays the dominant role for the deeply bound s states.!^ The s-wave optical 
potential for the symmetric nuclei pn = Pp is written as, 

2pVsir) = -47r[ei6op(r) + e2Bop\r)]. (2-7) 

For all 6o and Bq values satisfying the relation 

60 H PeBo = constant, (2'8) 

£1 

the optical potential ()2-7p has the same strength at p = pe and provide almost same 
structure of the pionic atoms. Actually the correlations between potential parameters 
are found phenomenologically and called as the Seki-Masutani (SM) correlations.'^ 
The SM correlations are expressed as; 

bo + asBo = (3s = i-Om + OMi)m-\ Os ~ 0.23m3, (2-9) 

and 

= /3p = (-0.2 + 0.02i)m-^ Op ~ 0.37m^, (2-10) 
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where, 

4 
7 = 1 + -vrAcoPe. (2-11) 

The parameters 6o, -Bq) cq, and Cq determine the strength of the pion- nucleus optical 
potential shown in Eq. ()A-4p . Seki and Masutani found that the series of potential 
parameter sets, which satisfy the relations (j2-9p and (j2-10p . reproduce the observed 
data by the X-ray spectroscopy reasonably well. Namely, it is very difficult to fix the 
unique set of potential parameters from the data taken by the X-ray experiments. 
In the modern x^ analyses of the data, one can expect the unique determination 
of the potential parameters and, actually, one may find literatures which report the 
unique determination. However, there still remain the strong SM correlations, and 
the x^-values in these analyses have the long deep valley structure along the SM 
relations of potential parameters. 

We can find from Eqs. (|2-8p and (|2-9p as. 



Pe = -as, (2-12) 

and the value of as determined by the analyses of atomic data indicates, 

1 
Pe ^ ^Po- (2-13) 

Hence, it means that the SM correlation found phenomenologically is the consequence 
of the fact that the nuclear density probed by various pionic atoms observed by X- 

ray is always close to -po. This feature is also found for some of the deeply bound 

pionic states ,'21' 'ESJ and thought to be one of the robust features of pionic states. In 
the early stage of the exploration of deeply bound pionic atoms, this robustness was 

an advantage to deduce the pion properties at pe — -po, however, to make a step 

further, we need to investigate the possibility to deduce the pion properties precisely 

not only at pe ~ —po but also at various nuclear densities. If we can determine {qq)p 

for various densities, it is expected to provide important experimental information 
for the studies of the chiral condensate beyond the linear density approximation. 

§3. Numerical Results and Discussions 

In this section, we investigate how we can deduce the (qq) value at finite density 
precisely from pionic atom data. As already mentioned, (qq) was determined in 
Ref. Ej), however, the recent theoretical work provide a different formula including 
the Z* factor as introduced in Section [2Tl In addition, the (qq) value was determined 
only at /) ~ 0-Qpo because of the limited sensitivity of the observed pionic atom data 
to nuclear density, which is related to Seki-Masutani correlations to the potential 
parameters as mention in Section 12.21 Thus, we study theoretically these points to 
find the proper way to deduce the {qq) value for various densities from the future 
observations. We also pay attention to the possible simultaneous observation of the 
non-yrast 2s state with the deepest Is statc'^^''^^' 
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3.1. (d,^He) spectra in experiments 

We first show the comparison of the theoretical calculations with the latest data 
of the pionic atom formation in Sn isotopes. In Fig. [H we show the calculated results 
together with experimental data reported in Ref. 7). The instrumental resolution is 
assumed to be AE = 394 keV FWHM (Full Width Half Maximum) as same as the 
data.!^"* We added the constant background to the calculated results and then scaled 
the calculated spectra to reproduce the strength of the resonance peak of the Is state 
formation in the experimental data for this resolution. We also show the calculated 
results using the same theoretical model with the improved resolution AE = 150 keV 
FWHM which is expected to be archived in the planned experiments.'^''^ We find 
that the theoretical calculation reproduce the experimental data reasonably well and 
provide the reliable interpretation of the spectra, which is an essential foundation of 
the research activity to deduce the pion properties from the ((i,^He) spectra. And 
as we can see from the figure, we conclude that the pionic 2s state can be observed 
as a peak structure together with the Is state in the experiments with the better 
resolution AE = 150 keV FWHM. This result helps to propose the new experimental 
activity and motivate us to proceed the theoretical analyses reported in this article. 

Simultaneous observation of the non-yrast 2s state with the ground Is state is of 
essential importance in the experimental viewpoints. We expect smaller experimental 
uncertainty in the energy difference between the two states while larger ambiguities 
are usually associated in the calibration of the absolute excitation energy. 

In order to resolve the two states Is and 2s in the excitation spectrum, the 
spectral energy resolution AE must be improved by a factor of about two from ~400 
keV to ~200 keV. The resolution is governed by two major contributions, namely 
the incident beam energy spread and the energy loss in the target. The latter can be 
controlled by changing the target thickness. However, adaptation of thinner target 
results in smaller statistical merits or in longer data accumulation periods unless the 
beam intensity is increased as compensation. 

A new experimental approach is thus started^ in the RI beam factory, RIKEN^^I 
to achieve improved resolution of 200 keV with much higher statistics. The experi- 
ment makes full use of the very high intensity deuteron beam of ~ 1 x 10^^ /second, 
which is more than 30 times higher than the previously available intensity in the SIS- 
18 accelerator in the GSI. Thus, the statistical merit is still huge even if the adopted 
target thickness of 5 mg/cm^ is factor of three smaller than the previous experiments, 
and the larger luminosity is essential for systematic study over wide range of nuclei. 
The other contribution in the energy resolution, the incident beam energy spread, is 
~ 5 times larger in the RIBF. Elaborate study on the accelerator, the beam transfer 
line, and the spectrometer is ongoing to establish a dispersion matching beam optics 
where the contribution of the beam energy spread is suppressed in the excitation 
spectra.'^ 

Simultaneous observation of the 2s state with the Is state is important also 
in the theoretical viewpoints. As shown in detail in Appendix [Bl the uncertainties 
of the theoretical calculation of the pionic atom formation cross sections due to 
the neutron wave function and the nuclear excited levels could be an obstacle to 
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Fig. 1. The (d,^He) spectra for the formation of the deeply bond pionic states in Sn isotopes as 
indicated in figure. The solid and dashed lines denote the theoretical calculations, which include 
only bound state contributions, while the histogram shows observed experimental data.l^ The 
instrumental energy resolution is assumed to be AE — 394 keV FWHM for the dashed lines and 
150 keV FWHM for the solid lines. The calculated results are scaled to reproduce the strength 
of the Is state formation peak in the observed spectra for AE — 394 keV FWHM. 



deduce the properties of pion by the precision measurements. Thus, as a practical 
way to reduce the uncertainties and deduce the rehable results, we should make use 
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Fig. 2. Calculated ^^^Sn(d,''He) spectra for the formation of the pionic bound states at Td — 
SOOMeV are shown as functions of the reaction Q-value, where the Fo factors in Table |V] are 
used. The dominant subcomponents are also shown as indicated in the Figure. Instrumental 
energy resolution is assumed to be 150 keV FWHM. The vertical line indicates the threshold 
Q = -141.6 MeV. 



of the recoilless kinematics to populate the plural number of pionic states coupled 
with the same neutron hole state as dominant contributions in the ((i,^He) spectra. 
And by comparing the strength of these contributions, we can effectively remove the 
ambiguities due to the structure of the target nucleus, the nuclear excited levels, 
and the neutron wave functions. In this context, the simultaneous observation of 
Is and 2s pionic states coupled with the (3si/2)^^ neutron hole is very interesting. 
We show the calculated ^^^Sn((i,^He) spectra in Fig. [2l This reaction is proposed in 
Refs. I14|) .ll5 p . As we can see form the figure, the Is and 2s pionic states coupled 
with the same (3si/2)n^ neutron hole state can be seen as the clear peak structures 
with the reahstic energy resolution AE = 150 keV FWHM. Thus, the ^^^Sn((i,^He) 
reaction can be one of the good reactions for our purpose. The expected spectra for 
Te isotopes are summarized in Appendix iDl 

3.2. Determination of bi parameter at various nuclear density 

As shown in Ref . JTj) , the observed binding energies and widths of bound pions 
have been used to determine the strength of the s-wave isovector potential param- 
eter (6i). The important points we should address here are to find out the way to 
determine 6i value at various densities to know its density dependence beyond the 
linear form. 

First, we consider the effective nuclear density probed by atomic pion.^^' The 
effective nuclear density pe is defined as the nuclear density at the radial coordinate 

Te as Pe 



p{re 



where the overlapping density, 
S{r) = p{r) \ R^l{r) 



|2 2 

r\ 



(3-1) 



has the maximum value. Here, Rni{r) is the radial wave function of the pionic atom 
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Fig. 3. Calculated peak positions r^ (upper frame) and the corresponding effective nuclear densi- 
ties pe (lower frame) for the observed pionic atom states by X-ray experiments are plotted as 
functions of nuclear mass number A. The quantum numbers of each atomic state are indicated 
in the figure. 



in a state of (n/). The definition of the overlapping density S{r) indicates that it is 
expected to evaluate the sensitivity of the iso-scalar s-wave optical potential term to 
the energy eigen values in the sense of the first order perturbation theory. In Fig. [3l 
we show the calculated r^ and pe systematically for the pionic atom states which have 
been observed by the experiments of the X-ray spectroscopy, so far.l^ As shown 
in the figure, r^ increased with the nuclear mass number A monotonically, however 
Pe is almost constant ('saturated') except for the very light nuclear cases. These 
features, which are almost independent on the quantum numbers of the states, can 
be understood by the pocket structure of the potential for the atomic pion and the 
localization of the overlapping densities there as can be seen in Appendix O Thus, 
all pionic atoms observed by the X-ray experiments probed the almost same nuclear 
density nearly independent on the nuclide and the quantum numbers of the pion 
except for the pionic states in the light nuclei. Thus, it seems that all observations 
of these states only provide pion properties around O.GpQ. To obtain pion properties 
at various nuclear densities, it will be necessary to observe other states. 

Then, we consider the atomic states in ^^^Sn up to n = 6 including the deeply 
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bound and non-yrast states, which can not be observed in X-ray experiments. When 
we use the Coulomb wave functions for bound pions, pe distributed in relatively 
wide range 0.05 ^ Pe ^ 0.16 [fm~'^] depending on the bound states. As naturally 
expected, the states with smaller orbital angular momentum i tend to probe the 
larger nuclear densities. As the realistic cases, we show the results calculated with 
the optical potential in Fig. [H pe and r^ only change inside the smaller range than 
those calculated with Coulomb potential only. This behavior can be also understood 
by the picture of the potential pocket and the localization of the overlapping densities 
at the nuclear surface. Thus, we find that the effective nuclear density is rather stable 
(0.08 < Pe ^ 0.11 [fm~^]) even for the non-yrast states and the deeply bound states. 

We mention here that an interesting tendency appeared in Fig. U] that the states 
with higher n for a fixed i provide larger pe values contrast to the usual intuitions. 
This is because of the larger tunneling effects to the central soft core for more lightly 
bound states. Thus, the sensitivity of pion in the higher n states moves to higher p, 
however at the same time, the absolute magnitude of the strong interaction effects 
are reduced rapidly as n increases. For example, we may think that it is better to 
observe the 4s and 4/ states to probe pion properties in different nuclear densities 
from Fig. HI This will be wrong since the strong interaction effects are too small for 
both states and the p-wave interaction plays dominant role in / states^ and hides 
the 6i effects even they may provide information for different p. 

Hence, we find that the observation of the binding energies of the pionic atoms 
will provide the pion properties near p ~ O.GpQ. Thus, to deduce the density depen- 
dence of the 6i parameter, we require the extremely higher energy resolution data 
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Fig. 5. Contour plots of the binding energies (left) and widths (right) for the Is (upper) and 2s 
(middle) states, and the difference of Is and 2s states (lower) in ^^^Sn in the &o-R.e(Bo) plane. 
The dashed lines depicted by SM are the parameter sets satisfying the Seki and Masutani 
correlation. ^ The numbers in the figure indicate the values of the binding energy and width 
of contour lines in unit of keV. 



as we will see later in the discussion with energy contour plots in this section. 

We, then, consider the contour plots of the binding energies B.E. and widths 
r of Is and 2s pionic states in vr"— ^^^Sn system. By the contour plots of eigen 
energies, we can see the difference of the sensitivities to the nuclear density as the 
deviations from the SM correlation and we can see the required energy resolution to 
distinguish the pion properties at different p. The contour plots of the differences 
of the binding energy and width of Is and 2s states are also shown, which could be 
used to deduce the systematic errors due to the calibration of the absolute binding 
energies and the uncertainties of the neutron distribution of ^^^Sn. We mention here 
that the binding energies and widths of the Is states in Sn isotopes were precisely 
determined in the last experiment. ' For example, the binding energy and width of 
the Is state in ^^SSn are B.E.= 3.744 ±0.018 [MeV] and F = 0.341 ±0.072 [MeY]^ 

We show the numerical results in Fig. [S] in bQ—HeBo plane together with the SM 
correlation line Eq. ()2-9p . We find that the contour lines for the binding energy of 
Is state are almost parallel to the SM line. Since the slope of the contour lines in 
bo-KeBQ plane provide the information of the nuclear density probed by the atomic 
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pion as shown in Eqs. ()2-9p and ()2-12p . this feature indicates that the nuclear density 
mainly probed by the Is state is very close to as as shown in Eq. (j2-12p . and that 
the precise measurements of the Is state provide the medium effects of pion at ag- 
in addition, it is very hard to determine the unique parameter set {bo, ReSo) on the 
SM line only from the observation of Is state. 

On the other hand, the contour lines of the 2s state show the slightly different 
slope from that of the SM line. This behavior indicates that the properties of the 
2s state are determined by the pion properties at slightly different nuclear density 
from as as indicated by the result shown in Fig HI And we may be able to find the 
unique parameter set (6o, Rei?o) by using the both B.E. of Is and 2s states with 
high precision. The plots of the difference of the binding energies of Is and 2s states, 
which is expected be observed with high precision without systematic errors due to 
the absolute energy calibration, show the similar behavior of that of the Is state 
because the larger binding energy value of the Is states dominate the behavior of 
the plot. We may also expect that the uncertainties due to the neutron distribution 
are partly cancelled in the plots of the energy differences. 
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In principle, to determine the 61 value for different p by observing the atomic 
states, we need to distinguish the nuclear density observed by these states, which will 
be equivalent to fixing the unique set of {60, Rei?o} by data, and then to determine 
61 value independently for each state. Thus, the experimental data with very high 
precision are necessary to obtain 61 values at various p as indicated in Fig. [5l 

The contour plots of the widths show different behavior from those of the binding 
energies, and the contour lines are not parallel to the SM line for all three cases shown 
in the right panels in Fig. [5j To understand this behavior intuitively is a little difficult 
since this is the effect to the imaginary eigenvalues from the modifications of the real 
part of the potential. However, we can naively expect that the p^ behavior of the 
imaginary potential can provide the different sensitivity of the widths of the pionic 
states from that of the binding energies. The contour plots of widths indicate that 
the precise determinations of the pionic widths will provide the constraints to the 
potential parameters, though the determinations of the widths are more difficult in 
general than those of the binding energies. 

In Fig. [6l we show another contour plot in 61— Imi?o plane. The 61 parameter is 
one of the most important parameters and has the very close relation to {qq), and 
the Imi^o parameter is the leading term to determine the width F in the s-states. 
The contour of the binding energy indicates the reasonable independence of B.E. on 
lmi?0) while the contour of the width indicates the complexity of the behavior of 
the width which depends strongly both on IuiBq and 61 . The contour of the binding 
energy shows a clear relation between the accuracy of the binding energy data and 
the hi parameter determination. 

One possible way to deduce the 61 values for different p is the parameter search 
using the most precise experimental data of the atomic s states, which are most 
sensitive to the s wave potential parameters ,'121 allowing the different value of 61 
for each states. Since this procedure, in principle, require to distinguish the nu- 
clear densities observed by the atomic states to determine the different 61 values for 
these states independently, the precisions of the data should be so high that we can 
clearly determine the 60 and Bq parameter uniquely and, thus, the feasibility of this 
procedure highly depends on the precision of the pionic atom data. 

3.3. (d^^He) spectra and observation of Z* 

The purpose of this section is to investigate the possibilities to determine the 
wave function renormalization factor Z* shown in Section [2. II from the observation 
of the formation cross section of the deeply bound pionic atoms. The basic idea is 
to observe the change of the cross section due to the modification of the pion wave 
function 4)1^ as, 

ct>iAr)^€ir) = Z:^"<t>iAr), (3-2) 

due to the renormalization factor Z* originated from the energy dependence of the 
pion selfenergy. The modification of the wave function causes the change of the the 
cross section as shown in Appendix \^ Z* is associated with the pion selfenergy 
in nuclear medium and, hence, has the p{r) dependence as shown in Eq. (|2-5p in 
contrast to that appeared in the standard text book of filed theory. Due to the p{r) 
dependence of Z*, we can expect to have different effects of Z* for the formation 
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Fig. 7. Pion radial density distributions |_R,r(r')p are shown as functions of the radial coordinate 
r for Is and 2s states in ^^^Sn. Solid and dotted lines show the pion density evaluated by the 
radial part of the renormalized {4>^) and unrenormalized (0) wave functions defined in Eq. ^2-2\ . 
respectively. The vertical line indicates the nuclear half radius R = 5.4761 fm. 



Table I. The calculated effective numbers of [(ls),r ® (3si/2)^^] and [(2s)^ ® (3si/2)n^] subcompo- 
nents by renormalized and unrenormalized pion wave functions defined in Eq. (|3.2p in PWIA. 
The ratios of the effective numbers of Is and 2s pionic states are also shown. 



pion wave function 



iVcff(ls) 
1.36 X 10~' 
1.41 X 10"^ 



jVeff(2s) jVeff(ls)/Afcff(2s 



2.92 X 10"^ 
3.00 X IQ-^ 



4.66 
4.70 



Table II. Same as Table U except for the results in DWIA. 



pion wave function TVeff (Is) 






iVeff(2s) iVeff(ls)/Afeff(2s) 



2.05 X 10" 
2.08 X 10' 



4.09 X 10' 
4.17 X 10" 



5.01 
4.99 



cross sections of different subcomponents [^tt ® Jn^] i^ general. 

We sliow in Fig. [7]the pion radial densities of Is and 2s states in ^^^Sn calculated 
by (j)^ and (j) in Eq. ()3-2p . Because Z* is one outside the nucleus, the pion densities 
are modified only inside the nucleus < r < 8 fm as can be seen in Fig \7\ In Fig. [51 
we show the contributions of (3si/2)n^ neutron-hole state to the ^^^Sn((i,^He) spectra 
for the formation of pionic atoms and we find the effects due to Z* factor are tiny. 
To estimate the effects of Z* to the observables, we show the calculated effective 
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Fig. 8. Contributions of the (3si/2)n^ neutron hole state to the ^^^Sn(d,^He) spectra for the for- 
mation of the pionic bound states are shown as functions of the reaction Q- value. The Fo and 
Fr factors are fixed to be 1. The solid line shows the calculated results with the renormalized 
pion wave function {4>^) and the dotted line that with the unrenormalized wave function (0) 
defined in Eq. (I3-2|I . The vertical line indicates the threshold Q = —141.6 MeV. 



numbers for the dominant subcomponents [(Is),^ ® (3si/2)ra^] and [(2s)7r ^ i'isi/2)n^] 
in Tables HI and HIl We find that the ratios Ncf[i^s)/Nes{2s) of N^d of pion Is and 
2s states formation, which are expected to be good quantities to deduce the pion 
properties independent on the uncertainties of the neutron wave function as discussed 
in Appendix [HI changes only around 1% for PWIA and 0.4 % for DWIA for results 
with (j) and </). These numbers seem to be too small to observe experimentally at 
present. Actually, the variation of the ratio of effective number N(,s{ls)/N(;f[{2s) due 
to the uncertainties of neutron wave function is larger and is around 10% as discussed 
in Appendix [Bj Hence, it seems difficult to deduce new information on Z* from the 
observed spectra by simply using the ratio of Is and 2s states formation strength. 
We mention here that the effects of Z* in the ratios of A'^efj are suppressed because 
the pionic Is and 2s states probe the almost same nuclear density as described in 
Section 13.21 and the Z* effects for both states are cancelled out in the ratio. The 
idea to extract new information on Z* from observables are considered to be still 
relevant. 

Finally, we also investigate the sensitivity of the neutron pick-up reactions to the 
optical potential parameters which satisfy the Seki-Masutani correlation Eq. ()2-9p . 
Since the effective numbers calculated by Eq. ()A-10p have different cpi^ dependence 
from the overlapping density Eq. ()3Tp . we may have chance to distinguish the 
potential parameters with Seki-Masutani correlation by the formation cross sec- 
tions. We consider three sets of 6o and Rc-Bq parameter which are, {bo, Rei?o} 
= {-0.0185,-0.05}, {-0.0300,0.00}, and {-0.0415,0.05}, where parameters are in 
pion mass units, and we find that there appear some discrepancies in the radial part 
of the integrand of Eq. (jATOp around < r < 4.5 fm. Since the nuclear half radius 
is taken to be 5.4761 fm in this case, the discrepancy only exists deep inside the nu- 
cleus which will be significantly suppressed by the distortion factor. Thus, we expect 
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that the effects to the formation rate are extremely smah. We have confirmed the 
expectation by numerical calculation which show that the variation of the ratios of 
the effective numbers of Is and 2s states for these potential parameter sets is less 
than 1 %. 

§4. Conclusion 

In this article, we have shown the newly calculated ((i,'^He) spectra on the ^^^Sn, 
^^"^Te, ^^^Te targets for the formation of the pionic atoms.'^''^ Based on these 
results, we have investigated the possibilities to deduce the new information on the 
{qq) value at various nuclear densities and on the wave function renormalization 
factor Z* paying attention to the recent theoretical and experimental developments 
of the studies of the pionic atoms. 

We have found that the formation spectra on the ^^^Sn target is suited for the 
observation of the pionic states because of the simple neutron level structure and 
large occupation probability of 3si/2 neutron state. However, the ((i,^He) reaction 
on Te isotopes could include the extra difficulties due to the complex neutron level 
structure and should be considered carefully for the experiments. 

We have also found that the nuclear density probed by the atomic pion is dis- 
tributed only inside the narrow region around p = 0.6/?o even for the deeply bound 
pionic states and the non-yrast states. This feature had an advantage in early stage 
of the exploration, however, it requires now the experimental data with excellent 
precision to deduce the information on {qq) at various nuclear density as shown in 
the contour plot studies in Section [3^ 

As for the uncertainties of the calculated cross sections and the determination 
of Z*, we have found that the ratio of the subcomponents coupled with the same 
neutron hole state is a good index relatively free from the systematic errors due to the 
neutron wave function. This observation is possible in experiments in RIBF/RIKEN. 
However, the renormalization factor Z* changes the pion wave function inside the 
nucleus slightly and the effects to the cross section are masked by the distortion 
factor. Since the nuclear densities probed by atomic pions are around 0.6/7o, the 
effects of Z* for various pionic atoms are almost same and are cancelled out in the 
ratio of the formation cross sections. Hence, the observation of Z* is rather difficult 
at preset. However, the idea to deduce Z* information from observables is important 
and the further studies are required. 

We think that we need to consider the pionic atoms in exotic nuclei to obtain 
the pion properties at various nuclear densities and the information on pion wave 
function renormalization. For example, the existence of the pionic nuclear states due 
to the strong Coulomb attraction^ and due to the thick neutron skin effects^ were 
predicted theoretically. Pionic atoms in unstable nuclei were studied in various cases 
in Refs. I35 ]) -I37 p . In these systems with exotic nuclei, we may have the different sen- 
sitivities of pions to nuclear densities and the different effects of the renormalization 
factor to the formation spectra. We will leave the study of these systems as future 
works. 
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Appendix A 

Theoretical Formula for Structure and Formation of Pionic Atoms 



The theoretical formula to calculate the structure and formation of the pionic 
atom3^™'EI'Tl3 are summarized in this Appendix. The energy spectra and wave 
functions of the pionic atoms can be obtained theoretically by solving the Klein- 
Gordon equation, 

[-V2 + ^? + 2//Kpt(r)] (t>{r) = [E- l^coui(r)]' ./-(r), (A-1) 

where // is the pion-nucleus reduced mass, E the eigen energy written as E = fi — 

Bj^ F with the binding energy B-,^ and the width F of the atomic states. The V^oui 

is the Coulomb potential with a finite nuclear charge density distribution Pchif)'- 



2 / Pch[r 



„'\ 



^coui(r) = -e' / p^dr'. (A.2) 

J \r — r \ 

The charge density distribution is written by the Woods-Saxon form as, 

Pch{r) = — r. "^^r, \/ T- (^■^) 

1 + exp[(r - Rch)/ach] 

The parameters of the charge distributions are taken from Ref. [38|1 and summarized 
in Table IIIII for the nuclei considered in this article. 

The V^pt in Eq. (jA-ip is the pion-nucleus optical potential, which we assume to 
be of the Ericson-Ericson type,'^ 

2fiVopt{r) = - 47r[6(r) + e2Bop\r)] + 47rV • [c{r) + e^^Cop''{r)]L{r)V, (A-4) 

with 

b{r) = ei[bop{r) + 6i[p„(r) - pp{r)]], (A-5) 

c{r) = ei^[cop{r) + ci[pn{r) - Pp{r)]], (A-6) 

L{r) = |l + ^7rA[c(r) + e.'CoP^r)]'^ , (A.7) 

where e-i and eo are defined as ei = l-\ and eo = 1 H with the nucleon mass 

M 2M 

M. As a standard parameter set, we use the potential parameters listed in Table HV] 
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Table 111. Radius parameters Rch of the charge distributions taken from Ref'38). The diffuseness 
parameters are fixed to be ach ~ t/4In3 for all nuclei with t — 2.30 fm. 



nuclide 


^i'^Sn 


^^"Sn 


^^^Sn 


i^^Sn 


i^^Tc 


l^^Tc 


Rch [fm] 


5.4173 


5.4588 


5.4761 


5.4907 


5.5368 


5.5617 



which are taken from Ref. [29j) . We use the Woods-Saxon form for the distributions 
of proton and neutron centers and assume the same shape for the both distributions, 

P{r) = PrXr) + Pn{r) = YT^^£-^J^^ (^-8) 

where R and a are the radius and diffuseness parameters, which are determined 
from the parameters Rch and ach of the charge distributions pch by the prescription 
described in Ref. H0|) . For the calculations of the bound states, we use the distribution 
parameters R and a same as the target nuclei of the ((i,^He) formation reaction. 
The pq is obtained by the correct normalization of the mass number of the daughter 
nucleus. 

Table IV. Pion-nucleus optical potential parameters'^ used in the present calculations. 



bo = -0.0283m;^' 


6i = -0.12m;;' 


Co = 0.223m~=' 


ci = 0.25m-=' 


Bo = 0.042im-* 


Co = Q.lOim-^ 


A = 1.0 





We can calculate the pionic atom formation cross sections in the effective number 
approaclP™ using the pionic atom wave function 0;^, the binding energy B.,^, and 
the width F obtained by solving the Klein-Gordon equation Eq. (|A-ip . The (d,^He) 
reaction cross section in the laboratory frame is expressed as. 



with 



^-ff = Y. \ [ drdax}{r)ei/2,mMMl ^ ^Ur,a)]jMXiir)f ■ (A-10) 

JMnis 

'^^^^' ( To ) indicates the elementary differential cross section at forward 

angles for the d-|-n — > ^He + vr" reaction in the laboratory system, which is extracted 
from the experimental data of the p -\- d ^^ vr^ -|- t reaction assuming charge sym- 
metry.EI'l3 AE is defined as AE = Q + m^ - B^ + Sn - 6.787 MeV for the {d,^Re) 
reaction with the pion mass m.,^, the pion binding energy B.,^, the neutron separation 
energy 5„, and the reaction Q-value. F denotes the width of the bound pionic state. 
For the neutron wave function -0^^, we adopt the harmonic-oscillator (HO) wave 
function in this article for simplicity. We also use the calculated neutron wave 
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function with Woods-Saxon type potential to check the theoretical uncertainties of 
the cross sections in Appendix [Bj We have used the oscillator parameter given by 
hw = 40A~3 MeV for the harmonic-oscillator wave functions, with A the nuclear 
mass number. The spin wave function is denoted as Ci/2,ms('^)) ^^d we take the spin 
average with respect to rus so as to take into account the possible spin direction of the 
neutrons in the target nucleus. Xi ^'^^ Xf express the initial and final distorted waves 
of the projectile and the ejectile, respectively. We use the Eikonal approximation 
and replace Xi &iid Xf according to 

x}ir)Xiir) = exp{iq ■ r)D{z, b), (A-11) 

where the distortion factor D{z, b) is defined as 

D{z,b) = exp - -cTdN dz' pA{z\b) - -GhN j dz pA-i{z' ,b) . {k-l2) 

Here, the deuteron-nucleon and ^He-nucleon total cross sections are denoted as (T^at 
and GhN- The function pA{z,b) and pA-i{z,b) are the density distributions of the 
target and daughter nuclei at beam-direction coordinate z with impact parame- 
ter b. The effective number approach is sometimes called the distorted wave im- 
pulse approximation (DWIA), and also called the plane wave impulse approximation 
(PWIA) in case we neglect the distortion effects in Eqs. (|A-lip and (jA-12p as putting 
D{z,b) = 1. 

In order to predict the spectrum of the ((i,^He) reactions, we need to take into 
account the realistic ground-state configurations of the target nuclei, the excitation 
energies, and the relative excitation strengths leading to the excited states of the 
daughter nuclei. To obtain a realistic total strength for the neutron pick-up from 
each orbital, we need to normalize the calculated effective numbers using the neutron 
occupation probabilities in the ground state of the target nucleus. The occupation 
probabilities are obtained from the analyses of the Z((i, t) Z reaction data and 
are not equal to one in general. 

As for the excited levels of the daughter nuclei, we use the experimental exci- 
tation energies and strengths obtained from the Z{d,t) Z reaction. Since the 
single-neutron pick-up reaction from a certain orbital in the target can couple to sev- 
eral excited states of the daughter nuclei, we need to distribute the effective numbers 
among these excited levels of the daughter nuclei in proportion to the experimen- 
tal strengths. Thus, the effective number for the pionic state {1-,^) formation with 
the A^-th daughter nucleus excited state coupled to a single neutron pick-up from a 
neutron orbit j„ is written 

A^cff(4 ® {3n'^)N) = N^s{i^0jn-^) X Fo{jn) X Fr{{j~^)n), (A-13) 

where A^eff(-^7r "^in"^) is the effective number defined in Eq. (jATOp . Fq the nor- 
malization factor due to the occupation probabilities of the neutron states jn in the 
target nucleus, and Fr is the relative strength of the A^-th excited states in the 
daughter nucleus coupled to the single neutron pick-up from the state jn- The Fq 
and Fpt of some medium heavy nuclei are compiled in Tables |Vl IVH \IX\ and |X] given 
in Appendixes [B] and \D\ and also in Refs.[TT|).ll2p. 
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Appendix B 

Nuclear Structure Dependence of Formation Cross Section of Pionic Atoms 



We discuss in this Appendix the uncertainties of the effective number approach 
used to calculate the cross sections. As we can see from the data in Ref. [Tj) and the 
theoretical prediction in Ref. I12p , the effective number approach works well to predict 
the shape of the ((i,^He) spectra, however, it fails to predict the absolute magnitude 
of the cross sections correctly. Thus, it is important to know the limitation of the 
applicability of this approach. This is also important to apply this approach to 
deduce Z* from experimental data as discussed in Section 13.31 

We show first the dependence of the calculated cross sections on the neutron 
wave functions V'n in targets. As mentioned in Appendix [AJ we adopted the HO 
wave functions for simplicity in this article. We calculate here the formation spectra 
using another set of ipn obtained by a theoretical potential Set OB in Ref. I41|) . 
The both neutron wave functions in 3,Si/2 state in ^^"Sn are shown in Fig. [9] as 
an example. The ^^"^Sn target nucleus was used in the latest experiment^ and 
considered in the theoretical calculation3l2i.lI2J before. The 3si/2 neutrons have 
dominant contributions to the cross section coupled with pionic s-states. We can 
see in Fig. [9] that the both wave functions show similar behavior, however the wave 
function calculated by the potential in Ref. 41) has a little longer tail than that of the 
HO. Because of the distortion effects, the long range tail part of the wave function 
can be important for evaluating the formation rate. 

Using these wave functions, we have calculated the ^^"^Sn (d,^He) spectra for the 
pionic atom formation and showed the results in Fig. [TOj The Fq and Fpi factors 
appeared in Eq. (|A-13p are taken from Table IV in Ref. [T2|) . We found that the cross 
sections calculated with ■0„ of Ref. HT]) is about factor 3 larger than those with the 
harmonic oscillator. The theoretical calculations in Refs. [TT]) . ll2p . which adopted 
the same harmonic oscillator wave functions for ipn with different parameterization 
of proton and neutron density distributions in K^pt, show the similar values for 
cross sections with the present case of HO wave functions, while the experimental 
result in Ref. |7|) shows the smaller strength than theoretical calculations. We think 
that the tail part of the wave function is significantly important to evaluate the 
formation cross sections because the inner part of the wave function is masked by 
the distortion effects. And at the same time, it will be very difficult to obtain the 
tail of ipn precisely and to make the accurate predictions of the absolute values of the 
formation cross sections. On the other hand, the shape of the spectra is relatively 
robust and is insensitive to il)^^ as shown in Fig. [lOl Thus, it is better to use the ratio 
of the formation rates of Is and 2s, for example, to deduce the information of Z* 
than to use the absolute peak height of each state to reduce the uncertainties due 
to tpn- The variation of the ratio of the effective numbers for the subcomponents 
[(ls)wC5i(3si/2)n^] ^-iid [(2s)7r<8)(3si/2)n^] ^^c to the different neutron wave functions 
is around 10 % for the case in Fig. [lOl 

We then consider the uncertainties due to the neutron occupation probability 
{Fq) in the target nucleus and the relative strength of the excited levels (-Fr) of the 
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Fig. 9. Neutron distributions of j„ — 3si/2 
state in ^^"Sn by the harmonic oscillator 
wave function (dashed line) and the calcu- 
lated wave function by the potential (set 
OB) in Refim (solid line). 



-140 -138 

Q [MeV] 

Fig. 10. Expected spectra of the ^^°Sn(d,^He) 
reaction for the formation of deeply bound 
pionic atoms at Td = 500 MeV plotted as 
functions of the reaction Q-value. The har- 
monic oscillator (dashed line) and the theo- 
retical (solid line) neutron wave functions^ 
are used. Instrumental energy resolution is 
assumed to be 150 keV FWHM. The verti- 
cal line indicates the threshold Q = —141.9 
MeV. 



daughter nucleus. We summarize Fq and Fr together with the excited energies (Ex) 
in Table |V] and IVD Fn and E^ of ^^^'^^SSn can be found in Table IV of Ref. USD. 
We show three sets of Fq values in Table |V] to estimate the uncertainties of the Fq 
values. Within these neutron states, it is known that the s^^ and d^/2 states have 
dominant contributions to the (d,^He) spectra in recoilless kinematics for the pionic 
atom formation. As we can see from Table |Vl the calculated results by RMF model 
with and without including the possibility of nuclear deformation (Fq and Fq) show 
the qualitatively same results. This fact indicates that these Sn isotopes are almost 
spherical nuclei. In addition, Fq obtained from experimental data also show similar 
results. On the other hand, the results of Te isotopes show significantly different 
and complicated features as shown in Appendix iDl Thus, the uncertainties due to 
Fq for Sn isotopes, which are roughly 20-30 %, are much smaller than those for Te 
isotopes. This ambiguity, however, does not affect the ratio between subcomponents 
with the same neutron hole state. 

As for the excited states of the daughter nuclei of the (d,^He) reaction, the 
relative strength Fr is listed in Table |Vl] for ^^iSn and in Table IV of Ref. |12]) for 
ii9,i23gj^^ It is difficult to evaluate the systematic errors of these numbers. We 
can say, however, that the strengths of the excited states of Sn isotopes are suited 
for the pionic atom formation because of their simplicity. The only 2^5/2 excited 
level of Ii9.i2i,i23gj^ isotopes has a little complex structure including a few levels, 
however, there is a clear dominant level with the largest value of Fr even for the 
d^/2- On the other hand the excited levels in I2i,i25rjig jj^clude several levels of the 
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Table V. Normalization factors, which correspond to the occupation probability of each neutron 
state for the ground state of the target nucleus, are shown for ^^"Sn, ^^^Sn, and '^^''Sn nuclei. 
The normalization factors indicated by Fo are evaluated by experimental datap2)H2 F^ and 
Fq are obtained by the theoretical calculations of RMF model!^ Fq is calculated by assuming 
the spherical shape of nuclei, and Fq by including the deformation effects. The numbers with 
asterisk for 137/2 and l/iii/2 states are estimated so as to satisfy the normalization of total 
neutron number in the valence shells by assuming the same occupation probabilities for both 
states. 



Fo Fg F^ 

Neutron ^^Ogn i^^Sn ^243^ Neutron ^^Ogn i^^Sn ^^^Sn Neutron ^^Ogn i^^Sn ^243^ 

orbit orbit orbit 



3Si/2 


0.70 


0.73 


0.80 


3Si/2 


0.54 


0.65 


0.78 


3Si/2 


0.53 


0.67 


0.78 


2d3/2 


0.50 


0.51 


0.69 


2d3/2 


0.64 


0.77 


0.85 


2d3/2 


0.68 


0.76 


0.85 


2d5/2 


0.94 


0.86 


0.94 


2d5/2 


0.93 


0.95 


0.96 


2d5/2 


0.94 


0.95 


0.96 


157/2 


0.70 


0.67* 


0.70* 


157/2 


0.95 


0.98 


0.97 


137/2 


0.97 


0.97 


0.97 


1/111/2 


0.42 


0.67* 


0.70* 


1/111/2 


0.28 


0.33 


0.47 


1/111/2 


0.24 


0.35 


0.45 



same quantum numbers with comparable strength as shown in Appendix [Dl which 
makes the reaction spectra rather comphcated and prevents us from deducing the 
information on pion. Thus, we can say that the Sn isotopes are more suited as the 
targets of the ((i,^He) reactions for the formations of the pionic atoms. 

Table VI. Excitation energy (Ex) and relative strength (Fr) of each excited level in ^^^Sn deter- 
mined from the experimental data of Ref. |4^|) . 



"^^T&r 



Neutron hole orbit E^ [MeV] Fr 

3si 



/2 
2^3/2 
24/2 



157/2 

lfell/2 



0.06 


1.00 


0.00 


1.00 


1.11 


0.65 


1.37 


0.35 


0.90 


1.00 


0.05 


1.00 



Appendix C 

Structure of Pionic Atoms in ^'^^Sn and the Effective Nuclear Density probed 

by Atomic Pion 

We show the calculated binding energies and widths in Table IVIII for vr^" atoms 
in ^^^Sn nucleus, which are expected to be observed in near future at RIKEN 
/RIBF.H^'ISJ We find, as the previous workP'l3'13''SJ showed, that the optical 
potential acts as the soft repulsive core and moved pionic wave functions outwards. 
This effect makes the widths of the bound states smaller to be quasi-stable states. 
As we can see from Table IVIII the eigen energies of pionic bound states are spread 
in wide region such as 250 < B.E. < 3800 [keV] and 10"^ < r < 320 [keV] for the 
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Table VII. Calculated binding energies B.E. and widths F oi tv — ^^^Sn atom in units of keV. 
B.E.FC indicates the binding energies calculated with a finite-size Coulomb potential only. 





^--i-Sn 


state 


B.E.FC [keV] 


B.E. [keV] 


r[keV] 


Is 


6227.7 


3829.2 


320.6 


2s 


I9I3.9 


I4I5.6 


76.8 


3s 


909.2 


732.5 


29.3 


4s 


528.3 


446.6 


I4.I 


5s 


344.7 


300.5 


7.8 


6s 


242.4 


215.8 


4.8 


2p 


2321.2 


2262.7 


II6.7 


3p 


1034.9 


I0I3.6 


39.6 


4p 


582.5 


572.9 


17.5 


5p 


372.8 


367.7 


9.2 


6p 


258.8 


255.8 


5.4 


3d 


1037.8 


1040.5 


2.5 


4d 


584.0 


585.6 


1.5 


5d 


373.6 


374.6 


9.0x10-1 


6d 


259.3 


259.9 


5.6x10-1 


4/ 


581.8 


581.8 


5.6x10-^ 


5/ 


372.5 


372.5 


4.8x10-^ 


6/ 


258.7 


258.7 


3.5x10-^ 



states with principal quantum number n = 1 ~ 6, and the shapes of pion densities 
are changed significantly for different states. 

We, then, consider the effective nuclear density probed by atomic pion, which 
was proposed in Ref. 28J and was one of the good quantities to know the nuclear 
density sensitively observed by pion in atomic states. The definition of the effective 
nuclear density pe and the overlapping density S{r) are given in Section 13.21 and 
Ref. [28|l . We show as an example calculated pion densities |i?P, nuclear density p, 
and the overlapping densities S for Is and 2s states in ^-^^Sn— vr" atoms in Fig. [TTl 
We find the same tendency as found in Ref. I28p for different systems. Namely, the 
peak positions of the overlapping density S{r) and, thus, the effective density pe 
values are almost same for Is and 2s atomic states even their eigen energies and 
wave functions are much different. We can interpret this tendency by considering 
the fact that the repulsive optical potential and the attractive Coulomb potential 
make the potential pocket at nuclear surface, where the overlapping densities are 
localized. 



Appendix D 

Numerical results for Te isotopes 



We show in this Appendix that the calculated results for Te isotopes, which are 
candidate nuclei of the future pionic atom experiments in RIBF/RIKEN.'i^''i3 In 
Table IVIIIl we show the calculated results of the binding energies and widths of the 
pionic states in ^^^Te and ^^^Te. We find that the level spacing is larger enough than 
the level widths and the states are quasi-stable as other deeply bound pionic atoms. 
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Fig. 11. Overlapping densities (lower frame) of the tt" densities (upper frame) with the nuclear 
density (middle frame) in pionic bound Is and 2s states in ^^^Sn. The vertical broken line 
shows the half-density radius of the nuclear density of ^^^Sn. 



In Table \1X\ we show the occupation probabihties (Fq) of the target nuclei 
^^^Te and ^^^Te. The factor Fq evaluated with the experiment is listed together 
with theoretical values, Fq and Fq. As you can see in the tables, the occupation 
of 3si /2 is significantly depend on the evaluation method. Since Fq is much smaller 
than Fq for the 3si/2 level, we can expect that both ^^^Te and ^■^^Te are largely 
deformed. However, the experimental Fq based on Ref . I44p are closer to the spherical 



nuclear value Fq. This feature seems somehow inconsistent, and these numbers could 
include large errors. Actually larger deformations for Te isotopes than Sn isotopes 
were reported in Ref. 45) . In Table (Xj we show the relative strength Fji for excited 
levels of daughter nuclei ^^^Te and ^^^Te. Here, we can see that the level structures 
of both nuclei are a little complicated. Especially, the 2^5/2 state splits into several 
levels, which include plural levels with similar strength. These structures may cause 
extra difficulties to deduce the pion properties from the formation spectra of the 
pionic atoms. 

We show in Fig.[T2]the calculated ((i,^He) spectra for the ^^^Te and ^^^Te targets. 
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Table VIII. Calculated binding energies B.E. and widths F of n 
units of keV. 



^Te and tt" 



Te atoms in 





-K — 


121^^ 


-K — 


i2bTg 


state 


B.E. [keV] 


r[keV] 


B.E. [keV] 


r[keV] 


Is 


4096.3 


372.5 


4037.1 


342.4 


2s 


1522.0 


93.3 


1507.3 


84.1 


3s 


789.1 


36.2 


783.4 


32.5 


4s 


481.6 


17.6 


478.9 


15.7 


5s 


324.2 


9.8 


322.7 


8.8 


6s 


233.0 


6.0 


232.1 


5.4 


2p 


2445.2 


141.9 


2432.7 


133.7 


3p 


1095.8 


48.4 


1091.4 


45.4 


4p 


619.4 


21.5 


617.5 


20.1 


5p 


397.6 


11.3 


396.6 


10.5 


6p 


276.7 


6.6 


276.1 


6.2 


3d 


1126.6 


3.3 


1126.7 


3.4 


Ad 


634.1 


2.0 


634.1 


2.1 


5d 


405.5 


1.2 


405.6 


1.2 


6d 


281.4 


7.5x10"^ 


281.4 


7.6x10"^ 


4/ 


629.5 


8.0x10"^ 


629.5 


8.5x10"^ 


5/ 


403.0 


6.9x10"^ 


403.1 


7.3x10"^ 


6/ 


279.9 


5.1x10"^ 


279.9 


5.4x10"^ 



Table IX. Normalization factors, which correspond to the occupation probability of each neutron 
state for the ground state of the nucleus, are shown for '^^^Te and ^^®Te nuclei. The normalization 
factors indicated by Fo are evaluated by experimental data.ES The normalization factors Fq 
and Fq axe obtained by the theoretical calculations of RMF model. ''^' The factor Fq is 



calculated by assuming the spherical shape of nuclei, 
effects. 



and Fq by including the deformation 



Fo 






n 






^0 






Neutron 


122 Te 


i26Te 


Neutron 


i22Te 


i26Te 


Neutron 


i22Te 


i26Te 


orbit 






orbit 






orbit 






3Sl/2 


0.34 


0.50 


3Si/2 


0.41 


0.69 


3Sl/2 


0.02 


0.05 


2^3/2 


0.31 


0.47 


2^3/2 


0.58 


0.80 


2^3/2 


0.36 


0.85 


2d5/2 


0.65 


1.00 


2d5/2 


0.90 


0.95 


2^5/2 


0.95 


0.98 


I37/2 


0.43 


0.59 


lff7/2 


0.96 


0.97 


lff7/2 


0.98 


0.99 


I/111/2 


0.23 


0.42 


I/I1I/2 


0.33 


0.51 


lftll/2 


0.41 


0.56 



We find that the shape of the calculated spectra strongly depends on the choice of Fq 
as we expected. The results with Fq and Fq show the qualitatively same behavior, 
however, the results with Fq show the much different behavior. This is due to the 
smallness of Fq for 3si/2 neutron state, which can have dominant contribution to 
the spectra coupled with the pionic s-states in the recoilless kinematics. The lack of 
the 3si /2 neutron contribution deformed the shape of the spectra drastically. In the 
spectra calculated with Fq, the contributions of many subcomponents compose the 
total spectra and, hence, they make it difficult to deduce pion information clearly 
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Table X. Excitation energy (Ex) and relative strength (Fn) of each excited level in ^^^Te and '^^^Te 
determined from the experimental data of Ref. I44[>. 



Neutron hole orbit E^lMeV] Fr 



3Sl/2 



2d; 

2dr 



3/2 

/2 



I57/2 
1^11/2 



0.00 


1.00 


0.21 


1.00 


0.48 


0.37 


0.59 


0.28 


0.92 


0.09 


1.17 


0.15 


1.32 


0.05 


1.36 


0.06 


0.45 


0.71 


1.15 


0.29 


0.29 


1.00 



i2brpg 

Neutron hole orbit E^^ [MeV] Fr 

3Si/2 



2d: 
2dr 



■3/2 

/2 



I37/2 
lfell/2 



0.00 


1.00 


0.04 


1.00 


0.64 


0.35 


1.05 


0.19 


1.14 


0.19 


1.27 


0.15 


1.43 


0.11 


0.64 


0.73 


1.74 


0.27 


0.00 


1.00 




-140 
Q [MeV] 



142 -140 

Q [MeV] 



Fig. 12. Calculated ^^^Te(d,^He) (left) and ^^''Te(d,^He) (right) spectra for the formation of the 
pionic bound states at Td = 500MeV are shown as functions of the reaction Q-value. Each line 
indicates the result with the different normalization factors {Fo,Fq,Fq) shown in Table HXl 
The vertical line indicated the threshold Q = -142.6 MeV (left) and -141.9 MeV (right). 



from the total spectra in this case. 

Thus, we can conclude here that the expected spectra of the ((i,^He) reaction 
on the Te isotope targets for the pionic atom formation mentioned in Ref. 14) could 
include large uncertainties due to nuclear structure which should be considered care- 
fully before experiments. 
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